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We gauge fix 600 SU(3) f3 — 6.0 configurations on a 16 4 lattice to a simple form of the maximal abelian gauge. 
We project the 5(7(3) valued links to the (7(1) x (7(1) subgroup, and extract (7(1) x (7(1) and monopole string 
tensions. After gauge fixing to the indirect center gauge, the (7(1) x (7(1) links are projected to Z(3) and a vortex 
string tension is measured. The vortex and magnetic current densities are measured. 



The idea that the long-range physics of an 
SU(N) gauge field can be compressed into the 
abelian sector of that field was suggested by 't 
Hooft Jl| . The simplest way to do this is to gauge 
fix the field so that it is as much in the abelian 
"direction" as possible and then to project to the 
abelian sector. 

The process of gauge fixing to the maximum 
abelian gauge (MAG) has been studied exten- 
sively in SU(2) lattice gauge theory. Gauge fixed 
configurations are projected to the (7(1) sub- 
group. From this (7(1) gauge field, monopole cur- 
rents are located and a potential from the (7(1) 
gauge field or from the monopoles can be found 
and compared to the full SU(2) case. Also stud- 
ied is the indirect center gauge (ICG). Here, after 
going to the MAG, the (7(1) field is also gauge 
fixed and then projected to the Z(2) center of 
the group There are still unresolved issues 
such as stability under cooling and sensitivity to 
gauge ambiguities for the case of SU{2) (3j. It is 
natural, however, to go on to the case of an 5(7(3) 
gauge group in order to look at confinement as it 
occurs in QCD. 

1. Maximal Abelian Gauge 

In SU (2) gauge theory in the continuum, the 
MAG is the gauge that minimizes the quantity 

C 2 ! - / [(4) 2 + (4) 2 ] 
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The equivalent lattice functional is 

C-l = £ Tr Ki^M*)^] . (2) 

A maximum of Eq. (|^) is reached by sweeping over 
the lattice and performing gauge transformations 
that maximize the lattice functional at a given 
site. 

In 5/7(3), the corresponding continuum func- 
tional is 

Gl u S = J [{Alf + {Alf + (Alf + (4) 2 + 

(Air + (Air] d* x . (3) 

On the lattice Eq.(|3|) becomes 

J2^[Ul(x)\ 8 U,(x)X s ]\. (4) 

For 5(7(3), gauge fixing updates are done based 
on SU(2) subgroups. We denote the three sub- 
groups as follows: SU(2) : Ai,A 2 ,A 3 , 5(7(2)' : 
A 4 ,A5,A 3 , SU(2)" : Xa,Xj,X'^, where A^ are the 
usual A-matrices, and A3 = diag( 1,0,-1) and 
A 3 ' = diag(0,l,-l). 

2. Abelian Projection 

Having reached a maximum of Eq. (|J) , it is nec- 
essary to project each 5(7(3) link matrix to the 
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£7(1) x £7(1) subgroup. In the SU(2) case, the 
projection to £7(1) is determined by the (oppo- 
site) phases of the diagonal elements of each link 
matrix. For SU(3), the phases of the diagonal 
elements of the link matrices do not determine 
a diagonal unitary 3x3 matrix. One simple 
way to project to £7(1) x £7(1) is to choose two 
of the diagonal phases, and determine the third 
by unitarity. Alternatively, each diagonal phase 
can be used, with unitarity enforced by an ad- 
ditional overall phase. The method used here is 
to find a diagonal matrix £7 such that the quan- 
tity Re (Tr \U A'~\ ) is maximized, where A is the 
SU(3) link. The purpose of this method is to find 
the "best" £7(1) x £7(1) approximation to a given 
5£7(3) matrix. Either of the two methods just 
described can be used for a trial solution. The 
trial solution is then successively updated using 
matrices of the form exp(j</>A3), or matrices with 
A3, A3 replacing A3. The procedure was found to 
give numerical accuracy in 10 iterations. As a 
side note, it would be the natural procedure to 
use for projecting spatially smeared £7(1) x £7(1) 
links back to the £7(1) x £7(1) manifold. 

3. Indirect Center Gauge 

With configurations projected to the £7(1) x 
£7(1) subgroup, it becomes a simple (and consid- 
erably less computationally intensive) matter to 
gauge fix to the indirect center gauge (ICG). This 
gauge is determined by maximizing the functional 

GIT = o^ElT^rf . ( 5 ) 
■Unlink X41 

where the £7 P are the C/(l) x £7(1) projected 
links. The easiest way to reach the ICG is with 
gauge transformations by factors of As, A 8 , and 
Ag, where A 8 = V3X 8 = diag(l, 1, -2), A 8 = 
diag(l, -2, 1), and A 8 ' = diag(-2, 1, 1). This 
choice of transformation has the advantage that 
solving for the optimal transformation gives an 
easy analytical solution. 

Projection to Z(3) is accomplished by choosing 
the element = 0,1,2, £= exp(27ri/3) that 

maximizes the quantity Tr^'lZ/t] = ^TrpTj], for 
each £7(1) x £7(1) link. This results in a set of 
Z(3) links. The Wilson loop is then calculated 



by multiplying these Z(3) valued links around the 
loop. 

4. Magnetic Currents 

Magnetic currents are extracted for each SU (3) 
color by applying the Toussaint DeGrand proce- 
dure to the U (1) x £7(1) links. This produces three 
magnetic currents, of which only two are indepen- 
dent. In our method j|, the sum over colors of 
the magnetic current is non-vanishing on a link 
by link basis. Nevertheless, when used to calcu- 
late a Wilson loop, this total current produces 
a null potential, as it should. After extracting 
the magnetic currents, monopole Wilson loops are 
calculated for each color. This part of the calcula- 
tion proceeds exactly as in £7(1) or SU(2) lattice 
gauge theory (|||. Finally the color average of 
these monopole Wilson loops is used to produce a 
monopole potential from which a monopole string 
tension is extracted. The color averaged fraction 
of links carrying magnetic current was also mea- 
sured and found to be 7.44(4) x 10~ 3 . 

5. Results 

We gauge fixed each of 600 configurations to 
the MAG using Eq. (|J) , generating one gauge copy 
per configuration. The MAG gauge fixing was 
done using over relaxation with overrelaxation 
parameter lu — 1.8. Gauge fixing was continued 
until <|X 12 | 2 + \X[ 3 \ 2 + ra 2 ) < 3 x 10- 11 was 
satisfied where the matrix X(x) is 

£ U^x)\ 3 Ul(x) + Ul(x - fi)X 3 U^x - n), (6) 

and the matrices X',X" are similar, but have 
A 3 , A3, respectively in place of A 3 . 

Further gauge-fixing using Eq.(j|) was per- 
formed to get to the ICG. The stopping criterion 
for this procedure was the same as in jjj, One 
gauge copy/configuration was generated, and the 
overrelaxation parameter for ICG was lu — 1.7. 
Links were then projected to Z(3) and Wilson 
loops calculated as in a Z(3) gauge theory. The 
P- vortex density, that is the fraction of Z(3) pla- 
quettes with non-zero flux, was measured to be 
2.25(1) x 10~ 2 . 
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Table 1 

Potential V(R) = aR + a/R+V a 



Potential 


a 


a 


V 


SU(3) 


0.053 


-0.267 


0.625 


U(l)xU(l) 


0.045(2) 


-0.07(2) 


0.19(1) 


ICG 


0.040(2) 


0.03(1) 


-0.02(1) 


Monopole 


0.037(1) 


0.02(1) 


-0.04(1) 



Wilson loops up to R — 8, T — 8 were gen- 
erated for the various cases. Potentials were 
determined from Wilson loop values by fitting 
- ln(W(H, T)) to a straight line in T. The ln(W) 
vs T plots were remarkably linear in all cases, 
even for T < R. From the potentials, string 
tensions were determined by fitting to the form 
V(R) = oR + a/R + Vq. Values for these pa- 
rameters are given in Table [l]. As can be seen 
there, the MAG [7(1) x C/(l), ICG vortex, and 
MAG monopole string tensions all lie below the 
full SU(3) value || by an amount well outside 
of error bars. In Fig.(|l|), we plot the potentials 
determined here along with the best SU(3) po- 
tential from the literature Q . The values of the 
constant Vq have been adjusted for display. 

Our findings for SU(3) are to be contrasted 
with the behavior in SU(2), where for one gauge 
copy/configuration, at couplings with compara- 
ble string tension, the MAG U(l) and ICG vor- 
tex string tensions are larger than the full SU(2) 
value H, while the monopole string tension lies 
very close to the full SU(2) value |(|. Account- 
ing for gauge copies with higher values of the 



V(R) vs. R for 1 6 4 ,5 = 6.0 



fuiiicliuuals causes all lliu string tensions Lu de- 
crease, with the £7(1) MAG sLmig leiisiun ending 

up ~ 8% below the full SU{2) result @. If the 
same trend holds here, the U(l) x U(l) MAG 
string tension would end up well below the full 
SU(3) value, with the monopole and ICG vor- 
tex values lower still. To remedy this situation 
for SU(3), we are studying more general forms of 
the MAG condition @ . 
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Figure 1. V(R) vs R for MAG (U(l) x U(l) and 
monopoles) and ICG (P- vortices) for SU (3) 
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